Shock Waves in Dispersive Eulerian Fluids

The long-time behavior of an initial step resulting in a dispersive shock wave (DSW) for the one-dimensional isentropic Euler equations regularized by generic, third-order dispersion is considered by use of Whitham averaging. Under modest assumptions, the jump conditions (DSW locus and speeds) for admissible, weak DSWs are characterized and found to depend only upon the sign of dispersion (convexity or concavity) and a general pressure law. Two mechanisms leading to the breakdown of this simple wave DSW theory for sufficiently large jumps are identified: a change in the sign of dispersion, leading to gradient catastrophe in the modulation equations, and the loss of genuine nonlinearity in the modulation equations. Large amplitude DSWs are constructed for several particular dispersive fluids with differing pressure laws modeled by the generalized nonlinear Schrödinger equation. These include superfluids (Bose–Einstein condensates and ultracold fermions) and “optical fluids.” Estimates of breaking times for smooth initial data and the long-time behavior of the shock tube problem are presented. Detailed numerical simulations compare favorably with the asymptotic results in the weak to moderate amplitude regimes. Deviations in the large amplitude regime are identified with breakdown of the simple wave DSW theory.